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We report on our preliminary results on the low-lying excited nucleon spectra which we obtain through a 
variational-basis formed with three different interpolators. 



There are two competing pictures for under- 
standing baryon physics in different mass regions. 
For hadrons with heavy quarks linear confinement 
and perturbative gluon exchange corrections P 
provide a satisfactory description. For lighter 
quarks towards the chiral limit we have to ex- 
pect manifestations of spontaneous chiral sym- 
metry breaking [2] . Lattice studies can shed light 
on the relevant mechanisms. In particular after 
the advent of Ginsparg- Wilson fermions it is now 
also possible to run computations relatively close 
to the chiral limit. 

However, exact chiral fermion actions (like the 
overlap action) are very costly. On the other 
hand, recently developed approximate solutions 
of the Ginsparg- Wilson equation allow to work 
at relatively small quark masses at only moder- 
ate numerical cost. Two such approximate so- 
lutions, the parametrized fixed point action and 
the cliirally improved action, have been studied in 
the BGR-collaboration [M4\ . Here we present re- 
sults for low lying nucleon states computed with 
the chirally improved (CI) action 5 . The CI op- 
erator is a parameterization of the Dirac opera- 
tor with terms essentially restricted to the hyper- 
cube. It allows for simulations at pseudoscalar- 
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mass to vector-mass ratios down to mps/my 
= 0.32 for a relatively large lattice spacing of 
a = 0.15 fm (see |3] for more details). Thus 
we can work with light pions on large physical 
volumes. The quark propagators are determined 
for quenched gauge configurations generated with 
the Liischer-Weisz action. Here we discuss results 
derived from 100 configurations on 16'^ x 32 and 
12^ X 24 lattices at a lattice spacing of a = 0.15 
fm as determined with the Sommer parameter [S] . 
This gives rise to spatial extensions of 2.4 fm and 
1.8 fm. In this study our quark mass values cover 
the region 0.02 < am^ < 0.2 corresponding to 
values of mps/my down to 0.397. A detailed 
account of our setting is given in [7] . 

We implemented the following three interpolat- 
ing operators for the nucleon 

Xi{x) = eabc[ul{x)Cj5db{x)]uc{x), (1) 
X2{x) = eabc[ul{x)Cdi,{x)]-f5Uc{x) , (2) 

Xsix) = i Eabc [uj (x) C '^Q-f5db{x)] Uc{x) . (3) 

These operators were Jacobi-smeared at both 
source and sink. We remark that we also ex- 
perimented with point sinks and wall sources but 
found unchanged results and no improvement of 
the quality of our data. We computed all cross 
correlations (n, m — 1,2, 3) 



(4) 



For a basis of infinitely many operators {n, m = 
l,2...oo) the diagonalization of the matrix C 
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would lead to the optimal operator combinations 
building the physical states. Finding these combi- 
nations is equivalent to the generalized eigenvalue 
problem ^ 

Cit)C^'\t) = AW(i,io)C(to)C^'HO, (5) 



with eigenvalues behaving as 



(6) 



Each eigenvalue corresponds to a different energy 
level Wk dominating its exponential decay. The 
optimal operators Xi which have maximal overlap 
with the physical states are linear combinations 
of the original operators Xi, 



E 



(7) 




0.20 0.25 
a m 

q 



Figure 1. Results for the nucleon masses as a 
function of the quark mass (all in lattice units). 
We compare data from 16'^ x 32 and 12'^ x 24 
lattices. 



The coefficients c^*'' are obtained from the j-th 
entry of the i-th eigenvector (corresponding to 
eigenvalue A'*-*). Note that here we work with 
only a finite basis of operators which gives rise 
to corrections to the single exponential decay Eq. 
© which are discussed in more detail in I7I8J . 

We project the operators in the correlation ma- 
trix to definite parity. This allows us to dis- 
entangle states with positive and negative parity. 
We use the eigenvalues of the full correlation ma- 
trix to identify the plateaus in the effective mass. 
In this region we then use fully correlated two- 
parameter fits for the largest and second largest 
eigenvalues to determine the energy of the ground 
state and first excited state. We also experi- 
mented with Bayesian fitting techniques but do 
not discuss these results here. 

In our data we clearly identify the nucleon, the 
two lowest negative parity states N(1535) and 
N(1650) and an excited state of positive parity 
with a mass well above the masses of the two ex- 
cited negative parity states and thus too high to 
be identified with the Roper state at 1440 MeV. 
We denote this state as N!^_. For negative par- 
ity and for the excited state in the positive par- 
ity sector the quality of the fits decreases towards 
smaller quark masses and we cannot maintain the 
standards of our fitting procedure. We do not give 
results in these cases. 



Our results for masses the are given in Fig. ^ 
Like previous studies |9llUllllI^ we clearly ob- 
serve a splitting between the nucleon and the neg- 
ative parity states. At large quark masses this is 
attributed to the orbital excitation of the quark 
motion in the confining color-electric field. The 
splitting slowly increases towards the chiral limit. 
This implies that near the chiral limit there is 
another mechanism contributing in addition to 
confinement. This is quite consistent with the 
chiral constituent quark model |2] where a consid- 
erable part of the splitting is related to the flavor- 
spin residual interaction between valence quarks. 
Also the splitting between the two negative parity 
states is clearly visible. 

At smaller quark masses (below a m,q — 0.05) 
plateaus in the effective mass plots for both neg- 
ative parity states and for N'_^ completely disap- 
pear and hence we cannot trace these states closer 
to the chiral limit. At small Euclidean time sepa- 
ration there is a hint from our data (for smeared 
source and point sink) that there may be one 
more positive parity excited state. However, our 
analysis did not provide consistent results for this 
weak signal. 

We performed our analysis for two different lat- 
tice sizes 16^ X 32 (full curves in Fig. ^ and 
12'^ X 24 (dashed curves). Fig. ^ shows that the 
two data sets fall on top of each other and we do 
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Figure 2. The mixing coefficients Ci of the op- 
timal operators plotted as a function of pseu- 
doscalar mass. By definition, the basis of the 
spectral decomposition has |cip + |c2p-l-|c3|^ = 1. 
Data from 16"^ x 32 is plotted with filled symbols, 
for 12^ X 24 we use crosses. 

not suffer from finite size effects. 

Interesting physical insight can be obtained 
from the mixing coefficients c^'-* ofEq. 0. In Fig. 
121 we show their dependence on the pseudoscalar 
mass and again compare IB'^ x 32 (filled symbols) 
and 12'^ x 24 (crosses) to illustrate the absence of 
finite volume effects. Operator xi clearly dom- 
inates the nucleon at all quark masses, whereas 
in QCD sum rule analyses 13 it is assumed that 
the equally weighted superposition of xi and X2 is 
the optimal combination coupling to the nucleon. 

The nucleon belongs to a 56-plet which would 
require that under the permutation of two quarks 
a diquark subsystem has positive parity since the 
color wave function is completely antisymmetric. 
In the large Nc limit the nucleon should not con- 
tain pseudoscalar or vector diquark subsystems 
and one expects the interpolator content from X2 
and X3 to be minimal as is indeed seen in our 
data. 

Fig. [3 shows that at large quark masses the 
A^(1535) is dominated by xi and iV(1650) by X2- 
As the chiral limit is approached 7V(1535) couples 
optimally to (xi — X2) while iV(1650) couples to 



(Xi +X2) and the X3 contribution is suppressed in 
both cases. iV(1535) and iV(1650) belong to the 
negative parity L = 1 70-plet of SU{6) and both 
contain scalar and pseudoscalar diquark compo- 
nents in their wave functions. So the mixing 
in the vicinity of the chiral limit is very differ- 
ent from the mixing in the heavy quark region. 
Within the quark model this mixing is attributed 
to the tensor force and our results hint at its re- 
lation to spontaneous chiral symmetry breaking. 

In 12 the subtraction of a parametrization of 
possible quenched artifacts (77'iV), giving rise to 
negative values of the correlation function for cen- 
tral i-values, was crucial for the identification of 
a Roper signal. We also observe such negative 
values but postpone a detailed analysis to future 
studies with better statistics. 
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